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Abstract 

Deviations between the form of trajectory assumed in a fit to a set of measurements 
and the actual form of the trajectory can give rise to sequential correlations in the 
residuals from the fit. These correlations can provide a more powerful goodness-of- 
fit test than that based on the minimum from a least squares fit. The use of 
this additional information is explored in the context of several common trajectory 
errors (e.g. decays in flight) encountered in charged particle tracking. 
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1 Introduction 



Trajectory fitting aims to determine a set of fitted parameters and to test 
the validity of the trajectory hypothesis. Both of these questions are usually 
addressed by minimizing a constructed from the squares of the deviations 
between the measurements and the parameterized trajectory. The value of 
at the minimum is used to test the adequacy of the fitted hypothesis. The 
test, however, explicitly ignores correlations amongst the residuals^ from 
the fit. Such correlations arise naturally in trajectory fitting, where adjacent 
measurements are in causal order. Deviations from the expected trajectory 

^ Email address: kowalews@uvic.ca 

^ The term residual in this paper refers to the signed distance of closest approach 
between a measurement and the fitted trajectory divided by the uncertainty assigned 
to the measurement. 
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due to either a discrete change at some point (e.g. a scattering or decay in 
flight) or to a continuous parameter change (e.g. dE/ dx energy loss or magnetic 
field anomalies) introduce correlated shifts in the positions of all subsequent 
measurements. The test is not very sensitive to these trajectory deviations 
when they are small on the scale of the measurement errors, since it considers 
only the squares of the residuals. 

This paper introduces new quantities for hypothesis testing that are applicable 
to any fits for which an ordering variable (e.g. time) can be identified.^ The 
mean correlation of ordered sets of residuals (nearest neighbor, next nearest 
neighbor, etc.) are used for this purpose. These correlations are essentially in- 
dependent of x^, and test the assumption that the measurements are mutually 
independent (or, more generally, that the correlations amongst the measure- 
ments are properly accounted for in the fit). It should be noted that the 
presumption of the independence of the measurements once known sources of 
correlation (e.g. multiple Coulomb scattering) are taken into account is also 
present in sequential fitting methods such as the Kalman Filter ^, and the 
correlation test developed here can be applied to the output of such a fit. 

The power of these correlation statistics for hypothesis testing is studied as a 
function of several trajectory deviations that arise naturally in charged particle 
tracking. For simplicity, the trajectories studied are circular arcs, correspond- 
ing to the projection of charged particle trajectories onto a plane transverse to 
a uniform axial magnetic field. The following sections describe the correlation 
variables, the simulation used to measure their effectiveness, and the improve- 
ment in discrimination between true and false hypotheses beyond what can 
be achieved using alone. 



2 Description of correlation variables 

The degree to which nearby measurements are correlated can be gauged by 
considering the mean correlation as a function of the distance between mea- 
surements. Using the standard correlation estimator to form an average 
correlation 



as a function of the distance between measurements gives fairly good discrim- 
ination between true and false trajectory hypotheses. The sums run from 1 

^ Not all fits are in this category: e.g. one cannot in general define a useful ordering 
variable for vertex fits. 
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to N — k, N is the number of measurements on the trajectory, Si is the signed 
distance to the fitted trajectory for measurement i, cXi is the estimated un- 
certainty of measurement i, and k is the correlation distance: k G [1,N — 1]. 
However, the following combination. 



Pk = 2^Wi / )_^Wi ; Wi = (2) 
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Fig. 1. The expectation value of the correlation variable pk is shown as a function 
of the correlation length k. The first 19 (pt) are shown for 20 (a), 40 (b) and 80 (c) 
measurements. The different symbols correspond to specific generated trajectories; 
the fitted trajectory is a circle in each case. The shaded region indicates the r.m.s. 
of the pk distribution for the circle trajectory, (d) shows (pk) for angular kinks 
occurring at different locations. 

motivated by considering the measure {6i — 5i+k)^ / {erf + crf_^^J, gives slightly 
better discrimination. This is because the correlation sought is in the actual 
distances S from the fitted trajectory, not in the residuals 6/ a. The weight fac- 
tor emphasizes those pairs of measurements with significant deviations. These 
correlation measures satisfy |pk| < 1 for all k. For uncorrelated measurements 
the expectation values of the pk are close to zero. Negative correlations are 
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introduced by the trajectory fit; they are small provided the number of fitted 
parameters is much smaller than the number of measurements. 

Fig. l(a-c) shows the expectation valueQ of each as a function of the cor- 
relation distance k (for k up to 19) for true circle trajectories and for three 
common trajectory deviations: a discrete angular kink, an uncorrected contin- 
uous energy loss and the decay in flight of a pion to a muon. The locations of 
the angular kink and decay in flight were uniformly distributed along the tra- 
jectory. The fitted trajectory in each case was a circular arc. The magnitudes 
of the (pk) for the incorrect hypotheses depend on the particular choice of 
parameters in the simulation. The number and location of the crossing points 
between (pk) for the incorrect trajectory hypotheses and (pk) for the correct 
hypothesis are very similar for the trajectory deviations studied here. Fig. 1(d) 
compares a sample with a discrete angular kink located near the center of the 
measurement region with a sample where the kink is located closer to the edge 
of the measurement region. The kink location clearly has an impact on the 
behavior of (pk), most notably for large k/N. 

Using all of the pk values gives the best discrimination between the correct 
hypothesis and a particular type of incorrect hypothesis, but does not pro- 
vide optimal discrimination for all types of incorrect hypothesis as can be 
seen by considering the different shapes (note in particular the locations of 
the extrema) of the trajectory hypotheses shown in Fig. 1. Furthermore, most 
of the discrimination power is concentrated at small k, since the r.m.s. of 
the Pk distributions expected for the correct hypothesis are smallest there. Q 
These considerations, along with the essentially linear behavior of the differ- 
ence (pk) false — (pk)true as a fuuctiou of the correlation distance k for small 
k/N, lead us to the following test statistic: 

A = E CkPk with Ck = 777^(L - k) and E = 1- (3) 

k=i ~ ^) k=i 

Choosing L as the nearest integer to N/8 was found to give good sensitivity 
to the trajectory deviations studied. 



3 Description of simulation 

The sensitivity of A to trajectory deviations was studied using a simple simu- 
lation. Charged particles were generated and tracked through a uniform axial 

^ Calculated numerically using the simulation described below. 

^ The number of correlation measures summed for a given k is N — k, giving a 

statistical error proportional to 1/ y/N — k. 
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magnetic field and measured points were generated in the plane orthogonal to 
the magnetic field direction. The measurement uncertainty was taken as Gaus- 
sian. Trajectories were generated with parameters typical of charged particle 
tracking detectors[||: 

• average hit resolution: 150 /im; 

• variation in resolution: factor of 5 between the best and worst measured 
points; 

• radial difference between first and last measurement layer: 54 cm; 

• number of measurements: varied from 20 to 160; 

• magnetic field: axial field of magnitude 1.5 Tesla; 

• initial particle momentum between 0.5 GeV/c and 5 GeV/c. 

The measurements were uniformly distributed along the trajectories, the hit 
efficiency was unity and no noise hits were generated. The location of dis- 
crete trajectory deviations (angular kink or particle decay) was randomly 
distributed along the trajectory. The generated points were fitted to a cir- 
cle assuming perfect pattern recognition, i.e. all generated measurements were 
used in the fit. 
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Fig. 2. The survival rate as a function of the location of a discrete angular kink 
within the measurement volume. The cuts on A and P(x^) are set to give 95% 
efficiency for true circular trajectories. The dashed lines bound the region used for 
subsequent study of the power of each variable. 
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4 Results 



The generated data were used to study the discrimination power of A and of 
as a function of specific trajectory deviations. Fig. 2 shows the fraction of 
trajectories surviving a cut on A or on the probabihty P(x^) as a function 
of the position of a discrete angular kink, where the kink angle was uniformly 
distributed between ±0.02 radians. The cuts were set to give 95% efficiency 
for true circular trajectories. There is little discrimination power for kinks 
occurring at either end of the measurement region. Based on this, only those 
trajectory deviations occurring within a fiducial region consisting of the central 
80% of the measurements were selected for subsequent study. 
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Fig. 3. The survival rate as a function of the size of a discrete angular kink within 
the fiducial region. The behavior for different N is shown in (a) for a cut that gives 
an efficiency e of 95% for true circular trajectories. The curves in (b) show the effect 
of varying the cut on e. 

The fraction of these selected trajectories surviving a cut that gives 95% ef- 
ficiency for true circular trajectories is shown as a function of the size of the 
angular deviation in Fig. 3. The correlation variable A is more powerful than 
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P(x^) provided N > 20 and becomes relatively more powerful as the number of 
measurements increases. This is to be expected, since the physical correlation 
length sampled by A is of the order of | of the track length, and increasing 
the density of measurements allows a more precise determination of the corre- 
lation. Similar behavior is observed for different choices for the efficiency for 
true circular trajectories. The power of A and P(x^) to discriminate vr — > 
decays from true circular trajectories is shown in Fig. 4. Note that the decay 
angle between the vr and fi does not in general lie in the measurement plane. 
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Fig. 4. The survival rate for vr — > decays occurring in the fiducial region as a 
function of the initial momentum of the pion in the plane transverse to the magnetic 
field. The behavior for different N is shown in (a) for a cut that gives an efficiency 
e of 95% for true circular trajectories. The curves in (b) show the effect of varying 
the cut on e. 

and that the momentum of the muon will be smaller than that of the pion. 
The correlation variable A is again a significantly better discriminant than is 

The extent to which the discrimination afforded by A is optimal was studied. 
The measurements along trajectories generated with a discrete angular kink 
were fitted using both the correct hypothesis, i.e. two circular arcs of con- 
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stant curvature with an angular deviation at a point, and using the nominal 
(incorrect) hypothesis of a single circular arc. The ratio of the resulting 
probabilities Pcircio(x^)/Pcorrect(x^) is an optimal discriminant variable in this 
case.[^ The position of the angular kink was taken as either the true position 
of the generated deviation (denoted "fixed Rkink" in Fig. 5) or as the best fit 
value after considering all potential kink positions (denoted "fitted Rkink 
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Fig. 5. The survival rate as a function of the size of a discrete angular kink for 40 (a) 
and 160 (b) measurements. The cuts on A and P(x^) are set to give 95% efficiency 
for true circular trajectories. Also shown are results from an optimal technique (see 
text) with the kink location either fixed to the generated position ("fixed Rkink") 
or determined by the fit ( "fitted Rkink" ) • 

Fig. 5). The kink position will in general be unknown, so the "fitted Rkink" 
curve is in practice the best one can achieve. As is seen in Fig. 5, the correla- 
tion variable A gives discrimination that approaches the optimal value as the 
number of measurements increases. 



° In practice, trajectory deviations of several different types may be present in a 
sample, so fitting for one particular type of deviation will not be optimal. 
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5 Discussion 



The correlation variables pk introduced here test the assumption that a set 
of measurements are mutually independent once known sources of correlation 
have been taken into account. Use of the simple combination of pk intro- 
duced here as A leads to an improvement in the sensitivity for detecting small 
trajectory deviations relative to that achievable using P(x^). This additional 
goodness-of-fit test is independent of the test, and can be applied to the 
residuals from both traditional least squares fits and from Kalman filter fits. 
The improved sensitivity to small deviations can be qualitatively understood 
by recognizing that has a quadratic dependence on individual deviations, 
while A has a linear dependence. For larger deviations (not shown here) P(x^) 
becomes a more powerful discriminant than A, as expected. For the deviations 
studied, the gain from combining the A and P(x^) tests was negligible. 

The A and tests have different dependencies on the input to the fit. The 
X^ test is sensitive to the scale of the assigned measurement errors cXi and 
can be compromised by mis-estimates of and non-Gaussian contributions to 
the measurement errors. The A test is insensitive to the scale of the It is, 
however, sensitive to correlations introduced in calibration procedures. The 
extent to which this is a practical problem in using the A test is a function 
of detector design and calibration. In particular we expect the A test to be 
most useful in devices where the effect of calibrations is randomized over the 
measurements on a trajectory (e.g. in small cell drift chambers, where the 
drift direction changes layer by layer), and to be less effective in devices where 
coherent effects dominate (e.g. in detectors employing a jet cell design). 

Potential uses for the A test in charged particle tracking involve selecting 
decays in flight and enabling high quality track samples (with reduced non- 
Gaussian tails on the track parameter resolutions due to trajectory deviations) 
to be selected. Given the ease with which it can be calculated, a test on A might 
serve as a filter for selecting tracks on which more computationally intensive 
tests [|, ^ will be performed. The correlations measured by the pk parameters 
may find use in a broader range of applications. 
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